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Abstract Onsager’s irreversible thermodynamics is used to perform a system¬ 
atic deduction of the kinetic equations governing the opening and collapse of 
transient pores in spherical vesicles. We show that the edge tension has to be 
determined from the initial stage of the pore relaxation and that in the final 
state the vesicle membrane is not completely relaxed, since the surface tension 
and the pressure difference are about 25% of its initial value. We also show 
that the pore life-time is controlled by the solution viscosity and its opening is 
driven by the solution leak-out and the surface tension drop. The final collapse 
is due to a non-linear interplay between the edge and the surface tensions 
together with the pressure difference. Also, we discuss the connection with 
previous models. 

Keywords Transient pore • Onsager irreversible thermodynamics • Spherical 
lipid vesicle • Membranes 


1 Introduction 

Biological pores are important entities in cells and play a crucial role in the 
spatial and temporal control of energy and matter fluxes. For instance, the 
modulation of membrane permeabilization mechanisms by the physicochemical 
properties of lipid bilayers is quite important for understanding antimicrobial 
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mechanisms m- When membrane organization is intact after transient pore 
formation, the permeabilization can be partially compensated by ion channels 
and pumps. When membrane organization becomes affected, this may eventually 
lead to the disruption of the membrane. In this case, the survival of the cell 
depends crucially of the speed with which the membrane is repaired 0. In 
general terms, pores in biological membranes have different origins and their 
mechanisms of formation may depend upon proteins and peptides E0- In the 
literature, it is suggested that transient pores are more commonly formed after 
the interaction between peptides and the lipid membranes m- Nevertheless, 
pore formation has also been studied in prebiotic environments, that is, without 
peptides and proteins, with the aim to understand their role as a transport 
mechanism in protocells modeled by unilamellar vesicles [5] . Exocytosis and 
endocytosis involve processes of membrane fusion and fission giving rise to the 
formation of a fusion pore , a channel, through which secretions are released 
from the vesicle to the cell exterior [M7MIMII]. The formation of the fusion 
pore occurs through the action of protein machinery (SNAREs) that have to 
apply mechanical forces in the range of 2-20pN for the opening of the pore 
[12]. Pores are also important for many technological applications such as 
nanomedice, sensing, and nanoelectronics US- 

Biological membranes are mainly composed of lipid molecules that form 
bilayers due to their amphiphilic structure. Lipid membrane are easily used 
in experiments to form uni- and multilamellar vesicles of different shapes and 
therefore they are very good models to investigate the main features of the 
formation and collapse of pores. In the simplest case, the pore formation and 
growth are generated by physical forces and do not require the presence of 
proteins. This process can be modeled by liposome rupture due to swelling. 

Recently, giant unilamellar vesicles immersed in an aqueous solution and 
composed of a single phospholipid bilayer and typical radius, R ~ 10 — 100/im, 
have been used to study the formation and collapse of transient pores pura 
1X81120) . The rupture and a transient pore appear as response to mechanics stress 
over the membrane. Once the mechanical stresses along the membrane exceed 
the lysis tension, the vesicle ruptures and a transient pore may appear [18] . 
Application of intense optical illumination on vesicles containing fluorescent 
membrane probes induces a sudden increase of the membrane surface tension 
until the creation of a pore in the bilayer Ematn]. Another technique used for 
creating pores in vesicles is the adhesion of the vesicle on an attractive surface 
m • More recently, Dimova and coworkers EMO] have developed a different 
method of inducing pores in giant unilamellar vesicles by the application of 
electric pulses. The dynamics of the pore is then captured by means of a fast 
imaging digital camera with a high temporal resolution. This method avoids 
some disadvantages of previous techniques, such as the use of fluorescent probes 
and glycerol to increase the solution viscosity, and allows better resolution of 
pore dynamics. For instance, the initial pore size is not well documented in all 
cases except in Ref. EBL for which the initial pore size is about r m i n = 0.75 pm. 
However, for vesicles with similar diameters, all these techniques produce pores 
with maximum radii in the range r max = 5 — 11 pm. 
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In this work we are interested in the theoretical description of pore dynamics 
after its formation. This dynamics can be divided into three stages. The first 
stage is characterized by rapid growth of the pore caused by the excess surface 
tension and the difference between the inner and the external pressures. The 
second stage is a relaxation process in which the slow shrinking of the pore is 
dominated by the edge tension. Finally, at long times the rapid pore collapse 
is due to a non-linear interplay among the main forces participating in the 
dynamics, the edge and surface tensions together with the pressure difference. 

The evolution of pores in membranes has been previously described in Ref. 
[2ll by adopting a model coming from the theory of viscous bare films in which 
pore opening is governed by the transfer of surface energy to viscous losses, and 
edge forces are neglected mm- An important assumption of the mentioned 
description is that the solution viscosity is small and the membrane viscosity 
controls the closing process. However, this seems to be in contradiction with the 
experimental fact that, in the case of membranes, the solution viscosity can be 
used to regulate the pore closing velocity [171118] . An extension of this theory 
was formulated in Ref. [23] where a lateral stress on the pore arising from the 
tangential movement of viscous aqueous solution relative to the membrane 
was included. This correction introduces a dependence on solution viscosity 
that accounts for the dependence of the pore closing velocity on the solution 
viscosity. As a consequence, the dependence of the pore collapse kinetics on 
the membrane viscosity becomes negligible in the initial stages. However, the 
final stage depends crucially on membrane viscosity. Other approaches taking 
into account the shape of the vesicle as a variable were proposed in [21] . 

In this article we use Onsager’s irreversible thermodynamics to deduce a 
model accounting for the kinetics of pore collapse in membranes. We begin by 
evaluating the free energy differential of an open vesicle, including the bending 
energy, the surface tension, the edge tension, and the coupling of the vesicle 
volume to its environment through a pressure difference term. Using the second 
law of thermodynamics we are able to deduce the kinetic equations that govern 
the time evolution of the pore and vesicle radii. By imposing consistent physical 
conditions on the resulting equations, we propose the correct dependence of 
the Onsager phenomenological coefficients on the pore and the vesicle radii. 
The obtained equations are numerically solved showing excellent comparison 
with experiments and a consistent behavior. In this way, we show that the 
closing velocity is mainly controlled by the edge tension and the solution 
viscosity. We also determine the time evolution of the surface tension and of 
the pressure difference. Our results question previous interpretations of the 
collapse dynamics. Previous models can be deduced from our general formalism 
depending on the election of the Onsager phenomenological coefficients. The 
corresponding inconsistencies are discussed. 

The article is organized as follows. Section [2] is devoted to discussing the ex¬ 
pression for the free energy differential of an open vesicle using thermodynamics 
arguments. Then we use the free energy differential and Onsager’s thermody¬ 
namics to derive the dynamical equations describing the process in Section [3] 
Finally, in Section |4j we discuss our results and present our conclusions. 
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Fig. 1 : a) Schematic representation of a vesicle with pore, b) Stressed vesicle and formation 
of a transient pore, c) Relaxed vesicle. 


2 Thermodynamics of a stressed vesicle 

Consider a relaxed vesicle radius R eq immersed in an aqueous solution at 
constant temperature T and total volume V = Vm + V ou t with Vi n = (4/3)7rf? 3 9 
the volume inside the vesicle and V out the volume outside the vesicle, see Fig. 
[T] c). The curved membrane has by definition a bending free energy per unit 
area given by Kb/R 2 q , with Kb = 2 k + R being k and R the bending and the 
saddle-splay moduli, respectively mm- Furthermore, due to the presence 
of the curved membrane, a finite pressure difference between the inside and 
outside of the vesicle exists, AP = P in — P out > 0. 

After an initial perturbation that increases the membrane surface tension 
to (Jo and the vesicle radius to i?o, a pore radius r{t) is produced that initiates 
a time dependent process in which the radius of the vesicle becomes a function 
of time, R(t), see Fig. [ljj. Here, we will analyze the dynamical process of the 
pore growth and collapse in terms of r(t) and R(t). 

The free energy change of the system has four contributions: dF = dFy + 
dF a + dFs + dFi, coming from volume (Fy), surface tension (F a ), bending 
curvature ( Fb ) and edge tension (F)). Thus, in the stressed state we have that 
the free energy has the following expression m 

dF = -APdV in + ddA + 7 dl , (1) 

where at the right hand side we used the relation dV ou t = ~dVi n and d 
stands for the total effective surface tension given by d = er + Kb/R 2 . The 
last term, responsible for pore closure is the contour free energy with 7 the 
edge tension and l = 2nr the pore contour length. In this form, for small pore 
radius ( r(t)/R(t ) « 1) the open vesicle volume and area may be written as 
Vi n = 47 rf ? 3 /3 and A = 47tR 2 — 7 rr 2 respectively, thus allowing us to express 
the free energy differential in terms of radii r and R: 

dF = 4 ttR (2 d — RAP) dR + 2n (7 — dr) dr. (2) 

Notice that before the initial perturbation that yields the appearance of a 
surface tension (er = 0 ) in the membrane of the vesicle, the radius of the vesicle 
is R eq . In that case, Eq. (J 2 J) reduces to: AP eq = 2Ks/Rlq- 

The explicit expression of the surface tension can be deduced by noticing 
that the energy required to stretch a membrane by reducing the thermal 
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fluctuations on the vesicle shape is given by the expression 


128] 


F a = O, 


{A — A eq ) 


2 A 


eq 


( 3 ) 


where o c is a characteristic surface tension, while A eq and A are vesicle areas 
before and after increasing the membrane tension. Since the temperature and 
the volume of the solution are controlled, the change in surface energy dF a 
with respect to the area of the membrane is 


a = 


(§) -*.(± 
\ J T,v \ 



( 4 ) 


Using the explicit expression of the area in terms of r and R the last expression 
takes the following form 


~ R 2 r 2 ' 

R\ q 4 R? eq ' 


( 5 ) 


Equation ([5]) allows one to calculate the evolution of the surface tension during 
the entire process of opening and collapse of the pore. A similar expression was 
used in m to express the area occupied by the lipid when the membrane is 
stretched. 

After the pore collapse, the final radius Rf of the vesicle may be, in general, 
slightly larger than R eq (see Ref. m and Figure 2 below) and therefore the 
final value of the pressure difference contains an additional contribution, with 
respect to AP eq , coming from a residual surface tension 


AP f = 


2dy 

Rf 


2a f 2k b 

hf + hf 


( 6 ) 


This equation states that the Laplace relation gives the final pressure with a 
correction term due to the bending energy of the vesicle. For a lipid bilayer, 
the value of the surface tension just before the opening of the pore takes the 
value oo ~ 10 -5 lV/m whereas the bending modulus is about Kb ~ 10~ 2O J. 
This implies that the bending modulus is k ~ 10 -20 J at 18 degrees Celsius 
[531 whereas the saddle-splay modulus is approximately given by R ~ —an 
with a ~ 1 or less [MU. If the vesicle size is smaller than a critical size, 
Req < {Kb/o e q)~ 1 / 2 ~ 10 _2 pm, then the correction term due to the bending 
energy of the vesicle cannot be ignored. In the case of giant unilamellar vesicles, 
the bending energy effect on the final pressure is small. 

The pressure difference AP decreases as the pore collapses and the radius 
of the vesicle decreases because a solvent outflow. Previous simulations of this 
system indicate that during the collapse the vesicle shape remains approximately 
spherical [ 55 ] . Hence, an expression for the pressure difference can be obtained 
by calculating the mass outflow per unit time for a sphere with a pore which 
can be expressed in the general form 



pv • ds. 


( 7 ) 
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where p is the solvent mass density, v is the flux velocity crossing the pore with 
differential surface area ds. Using the Navier-Stokes equation for stationary 
states and incompressible fluids, rj s X 2 v — VP = 0, the approximate expression 
for the velocity in this case is 


Considering cylindrical coordinates and combining equations 0 and (|8|, we 
obtain the following formula for the pressure difference in terms of the mass 
outflow per unit time 



( 9 ) 


An explicit relation in terms of the vesicle radius will be given in the next 
section. 


3 Onsager’s irreversible thermodynamics and the dynamic 
equations 

The equations governing the dynamics for r(t) and R(t) can be deduced fol¬ 
lowing Onsager’s irreversible thermodynamics m that establishes phenomeno¬ 
logical relations for the time evolution of the system thermodynamic variables 
in terms of the conjugated generalized forces [27j. This procedure is systematic 
because it uses the second law of the thermodynamics through the calculation 
of the entropy production. 

From thermodynamics, it can be shown the general relation: dF = —TdiS < 
0, where d,S is the entropy produced during an irreversible process occurring 
at constant volume and temperature m- Therefore, the entropy produced 
per unit time is proportional to the time change of the system Helmholtz free 
energy 


diS dF dF dr dF dR 

dt dt dr dt dR dt 


( 10 ) 


In this expression, the time derivatives play the role of flows, Jk = d(k/dt, 
whereas the partial derivatives play the role of generalized forces Xk = — dF/d£k, 
m- Following the rules of non-equilibrium thermodynamics, we may assume 
that the flows are linear functions of the forces Jk = Yjj LkjXj, where the 
Lkj are the so-called Onsager’s phenomenological coefficients. Additionally, 
the elements Lkj also obey the Onsager reciprocal relations Lkj = Ljk . These 
reciprocal relations are associated with cross effects [32]. Therefore the linear 
laws imply 


at or on 




dt RR dR Rr dr' 


( 11 ) 
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The Onsager’s coefficients are mobilities characterizing the change of the 
radii in time (flows) to the change of the energy as a function of the radii (forces). 
Therefore, the units of the Onsager coefficients are the inverse of Jules-second 
per square meters: [L / 1 ] = J s/m 2 . Furthermore, the experimental results 
uzi suggest that the duration of the pore dynamics depends on the viscosity 
of the medium in such a way that the dynamics is slower when the solvent 
viscosity increases. Since the dimensions of the viscosity are [ 77 ^] = J • s/m 3 , 
then we have that 


L rr risk, ^rR — ^Rr Vsk and oc 77 ^ 3 , ( 12 ) 


with li (i = 1,2,3) a characteristic length of the system. In general, the 
Onsager coefficients may be functions of the state variables r or R [27]. The 
election of this dependence is crucial and a proper choice leads to a physically 
consistent description of the process. In the following we will show that pre¬ 
existing models can be obtained from our general theory by assuming constant 
coefficients, whereas a fully consistent system of equations need to consider 
that the coefficients should depend on the state variables. 

The change with time of the free energy given by Eq. © is therefore 

r] JT' d T? c\t 

—— = 47 ri? (2d — RAP) -j- + 2 - 7 T (7 — or) —. (13) 

at at at 

In this form, following the general scheme already presented we may write 
down the following set of dynamical equations for the time evolution of the 
pore and vesicle radii 


— = 2 - 7 T L rr {a 1 — 7 } — 47 t RL r R {2a — RAP }, 
at 


dR 

dt 


—AitRLpir {2d — RAP} + 2 'kLr t . {dr — 7 } . 


(14) 

(15) 


Equations (14) and (15) with a non-zero initial condition for the pore radius 
predict the dynamics of the transient pore in spherical vesicles. Both equations 
are coupled through the effective surface tension and the pressure difference 
terms. 


3.1 Dynamic equations and Onsager coefficients 


Dimensional analysis establishes that the Onsager coefficients have to be 
inversely proportional to a viscosity and a characteristic length. Irreversible 
thermodynamics states that phenomenological coefficients can in general be 
functions of the state variables. Therefore, here we will propose an expression 
for the Onsager coefficients as functions of the system’s state variables r and 
R, that yields a physically consistent model of pore dynamics. 

First, we have to emphasize that Eqs. (14) and (15) are not physically 


consistent if the Onsager coefficients are constants, because in this case when 
the pore collapses (r = 0) the time derivatives dr/dt and dR/dt are proportional 
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Fig. 2: Comparison of the numerical solutions (solid lines) of Eqs. and \l81 ), along 
with experimental data (symbols) taken from Ref. m-- (a) pore radius as a function of time, 
(b) vesicle radius as a function of time, (c) surface tension and (d) pressure difference. The 
parameters taken from the experimental measurements are R eq = 19.6pm, Ro = 20.5 pm, 
ro = 1.10 pm, cr c = 2.9 x 10 — 5 N-m~ 1 , r] s = 32 cP, and = 29 x 10“ 20 J. The edge tension 
was considered as a fitting parameter and was found to be equal to 7 = 0.92 pN. 


to — 7 , which has been considered in the literature as a constant parameter pn 
mrnmm- Hence, at r = 0 the dynamics does not vanishes. 

The Onsager coefficients are the inverse relaxation times of the pore dynam¬ 
ics. As implemented in experiments |16 |1171IT8] . changing the solution viscosity 
r] s allows to slow down the pore dynamics. Thus, we may assume that the 
Onsager coefficients scale with the solution viscosity r] s . In view of this, a 
consistent choice of the coefficients that satisfies the conditions ( p~2| ), does not 
introduces unknown free parameters and satisfies conditions dr/dt = 0 and 
dR/dt = 0 at r = 0, is given by the relations 


L rr — Lrr — 


1 


27 rr 0 r] s R 


and 


L r R — 


2nr 0 rj s R 2 ' 


(16) 


where 27rro is the initial length of the pore (perimeter) in which the initial 
radius ro can be estimated from Eq. ([5]). Note that both Onsager coefficients 
are well defined for all values of the radii. At the initial condition they take 
the values: L rr = Lrr = l/(2nr] s Ro) and L t r = 7 'o/( 27 r 77 s f?g), whereas for the 
final state they vanish. A more general expression for the Onsager coefficients 
can be proposed by considering the effect of the lipid viscosity 772 in the form: 
La ~ ( 77 s r 0 + 772 d ) _1 (r/f?), Li, ~ (r] s r 0 + 772 d ) -1 ( r/R) 2 with d the thickness 
of the bilayer. This second term could be of importance in the case of low 
viscous environments, see Ref. pg. 


Thus, Eqs. (14) and (15) take the form 
dr 

roVs-R = 


r or] s 


' dt 
dR 
dt 


: ( dr ~ ^ - 2 {2v - R A p } ^ 

(17) 

7* 2 

—2 {2<j — RAP} r + {dr — 7 } — 

R z 

(18) 


This set of equations must be complemented with the expressions for the 
effective surface tension d = a + Hb/R 2 with a given by Eq. ([ 5 ]), and the 
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Fig. 3: Comparison of the numerical solutions (solid lines) of Eqs. and \l8 1 ), along 
with experimental data (symbols) taken from US-' (a) pore radius as a function of time, (b) 
vesicle radius as a function of time, (c) surface tension and (d) pressure difference. The 
parameters taken from the experimental measurements are R eq = 19.7 pm, Ro = 20 Apm, 
ro = 1.15/xm, cr c = 3.0 x 10~ 5 N-m -1 , rj s = 1.13 cP, and = 29 x 10 “ 20 J. The edge tension 
was considered as a fitting parameter and was found to be equal to 7 = 1.21 pN. 


pressure difference AP given by Eq. |9|. The pressure difference AP is related 
with the time variation of the vesicle volume in the form 


Q = = -4:TrpR 2 R. 

at 


From equations ([9]) and (19) we obtain 


AP—***R. 


(19) 


( 20 ) 


Equations and ( |18| ) constitute a coupled set of two ordinary non-linear 
differential equations that have to be solved numerically along with ([5]) and 
(20). The results of the numerical integration of the radii r(t) and R(t) (solid 
lines) are shown in Figures [2|4] together with experimental data (symbols) 
taken from mmm- Our model reproduces very well the three stages of the 
dynamics of the pore radius (opening, decay and collapse). In addition, our 
model allows determining the time dependence of the effective surface tension 
and of the pressure difference as it is also shown in Figures [2]4j 

The values of the solvent viscosity, r] s , the equilibrium vesicle radius, R eq , the 
initial pore radius, r 0 and the initial vesicle radius, Ro indicated in Table [T] and 
used in Figures [2j [3] and [4j were taken from experimental reports from Refs. El, 
ID, and [2T], respectively. According to Ref. JM], we assumed that the bending 
modulus is equal to Kb = 29 x 10“ 2O J in all cases. The characteristic surface 
tension was estimated by using the expression cr c = ( 487 r k 2 ) / (R 2 q K bT ), 121, 
[24] . The last column in Table [l] reports the values of the edge tension obtained 
from the fits. 

In order to conclude this section, we want to discuss the assumptions 
underlying a very well know model which has been widely used for estimating 
the value of the edge tension 7 from experiments [21j . This model is based on 
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Fig. 4: Comparison of the numerical solutions (solid lines) of Eqs. 0 and \l8 along 
with experimental data (symbols) taken from Ref. m i (a) pore radius as a function of time, 
(b) vesicle radius as a function of time, (c) surface tension and (d) pressure difference. The 
parameters taken from the experimental measurements are R eq = 14.5 pm, Rq = 14.7 pm, 
7*0 = 1.15 pm, cr c = 3.1 x 10 — 5 N-m~ 1 , rj s = 32 cP, and = 29 x 10 —20 J. The edge tension 
was considered as a fitting parameter and was found to be equal to 7 = 0.57 pN. 

Table 1: Physical quantities used in the fittings of the experimental data reported in 


i)s(cP) Req(lim) ro(jum) Ro{/J.m) K b (10 20 J) a c ( 10 s N-m) 7 (pN) 


32 

19.6 

1.10 

20.5 

29 

2.9 

0.92 

1.13 

19.7 

1.15 

20.4 

29 

3.0 

1.21 

32 

14.5 

1.15 

14.7 

29 

3.1 

0.57 


the following equation for the pore radius 


dv 

— = 2TiL rr {err - 7 } 

and on the following equation for the vesicle radius 

3?7 S R 3 ' 


( 21 ) 


( 22 ) 


In writing this system of equations one has to assume that the vesicle free 
energy does not depends on the vesicle radius R and that the total surface 
energy only depends on the surface tension a [given by Eq. <§] , that is, the 
bending contribution Kb is neglected. In addition, the equilibrium Laplace 
relation between the surface tension and the pressure difference is assumed 
to hold during the whole process, AP = 2a/R. As a consequence of these 
assumptions, the model implicitly assumes that = 0 (no cross effects) 
and that L~^ = 4^772 d, with 772 the membrane viscosity (which is an unknown 

reproduce the 


parameter) and d the bilayer width m- Eqs. ( |21[ ) and 
experimental results provided that the vesicle radius dynamics, dR/dt , is 
controlled by the solution viscosity r/ s whereas the membrane viscosity, 772 , 
dictates the pore dynamics. This fact seems to be inconsistent with experimental 
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evidence (men. In addition, Eq. © does not satisfies the condition that it 
must be a critical point, that is, dr/dt = 0 and dR/dt = 0 at r = 0 and R = 0. 
From this equation it is clear that dr/dt | r=0 A 0 yielding to negative the pore 
radius. 


4 Discussion and conclusions 


In this work, we propose a theoretical model that accounts for all the features 
of the dynamics of transient pores in spherical lipid vesicles, see Figures]^ a), 
[Sja) and[4|a). Essentially, it allows to determine the value of pore edge tension 
along with the time course of the vesicle radius, the effective surface tension 
and the pressure difference from the measurements of the pore radius m 
fig] , see Figures [2jb)-(d),[3]Jb)-(d) and[4jb)-(d). All the parameters included 
have a clear physical meaning and can be directly measured in experiments; 
for example, the viscosity of the solution, the bending modulus or the initial 
surface tension, see Table |Tj 

The model was deduced by using Onsager’s irreversible thermodynamics 
and it consists on two coupled first-order time differential equations for the pore 
and the vesicle radii, r(t) and R(t) respectively. The coupling emerges from 
the dependence of the surface tension and the pressure difference on the state 
variables (r and R) and from cross effects coming from the linear laws of non¬ 
equilibrium thermodynamics. The latter introduce Onsager’s phenomenological 
coefficients whose dependence on r and R is provided to have a consistent 
physical model. This condition is not fulfilled by previous models [2Tll24j . In 
agreement with experiments, the characteristic relaxation time of the pore 
dynamics is controlled by the solution viscosity, T re i ax ccr) s , see Eqs. © and 


(18). 


The continuous lines in[2jb), [3^b) and[4jb) show the evolution of the radius 
of the vesicle, as predicted by the numerical solution of Eqs. |l7| and (18). This 


behavior is in agreement with the experimental results, Fig. |2[b). In particular, 
in experiments with constant illumination, this evolution is typical, and the 
vesicle radius may be reduced by 40% [IS]. The panels (c) and (d) of Figures [ 2 ] 
[3] and [4] show that after the pore formation, the surface tension and the pressure 
difference decrease drastically until the pore radius reaches its maximum value. 
After that, both increase slowly until a final equilibrium value of about 25% of 
their initial value. 

The previous results show that the membrane is not fully relaxed after the 
pore collapse. Furthermore, correlating the results shown in Figures [2|yi), [d|a ) 
and[4]ja) with their counterparts in (c) and (d), and with Eqs. ( |l7| and ( |18[ ), 
we see that it is more convenient to calculate the edge tension, 7 , from the 
slope of the initial quasi-linear relaxation regime, because the change of the 
slope in the final collapse stage manifests clearly the non-linear nature of the 
process attributable to the interplay among the edge and the surface tensions 
together with the pressure difference. In Fig. [2]ja) we show this through the 
average slopes of the two regimes, shown by the solid and dashed straight lines. 
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Near the initial quasi-linear relaxation regime, both a and AP nearly vanish 
[see Figs. 2(a)- 4(a)] and, therefore, Eqs. and ( fl8| may be approximated 
to give the following relaxation equation for the pore radius 


dr 

dt 


roVsRo 


(23) 


From this equation we can obtain a simple and useful expression for the edge 
tension as a function of the pore radius 


7 - -roVsRo 


A In I 


At 


(24) 


Using this expression we can calculate the values of the edge tensions reported 
in the first column of Table 2. These results correlate well with those predicted 
by the whole model that incorporates all the mechanisms participating in the 
relaxation dynamics; that is, it includes the effects due to the surface tension 


and the pressure difference. If Eq. (24) is used in the final regime (a rough 
approximation) then the values obtained for the edge tension are about one 
order of magnitude larger than in the first regime. Notice, in addition, that 
these values almost duplicate the corresponding ones reported in previous 
works mmm, see the third and fourth columns of Table 2. The reason is 
that previous estimations of the edge tension have been obtained using the 
model proposed in Ref. |2Tj. As discussed after Eqs. (21) and (221 , the election 
of the Onsager coefficients in that case corresponds to L~} ~ Airi^d whereas 
in our case should be L~^ ~ 27rro7 s , involving a factor of two. 


Table 2: Comparison of the edge tensions 7 obtained in the different relaxation regimes 
and with previous results reported. 


Figure 

Initial regime 

Whole model 

Final regime 

Previous results 


2 

3 

4 

1.60 pN 

1.32 pN 

1.01 pN 

0.92 pN 

1.21 pN 

0.57 pN 

26.5 pN 

14.3 pN 

10.8 pN 

10.0 pN, Ref. 02] 
14.3 pN, Ref. 01] 
6.0 pN, Ref. [21] 


The values of the parameters taken from experimental data reported in 
the literature are the solvent viscosity, the initial and final vesicle radius, the 
characteristic value of the surface tension and the bending modulus mm 
mmm- The magnitude of the initial pore radius can be obtained from 
the relation for the surface tension. Therefore, all fits were obtained by only 
adjusting the edge tension. 

It is convenient to mention that pore formation in vesicles immersed in high 
viscosity solutions (32cP), Figs. 2 and 4, is about Stf orrnat i on = 0.1 —0.2 seconds. 
In contrast, when the solution’s viscosity is near that of the cytoplasm (bewteen 
1-3 cP, as in Figure 3), then the formation time is much less, 6tf orrna tion = 
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0.0075 seconds. This second formation time scale is in the order of magnitude 
with what is expected for pore dynamics in living cells. It is interesting to 
notice that for the low viscosity experiment, Figure 3, all the values reported 
in Table 2 are consistent between them, in contrast with the data obtained in 
the other two cases. 

To conclude, we may stress that Onsager’s irreversible thermodynamics is a 
powerful formalism that allows to deduce simple and consistent physicochemical 
models for biological small systems, in particular, to describe the dynamics 
of pores in biological membranes. The obtained equations are consistent with 
the second law of the second law of Thermodynamics and it may easily be 
generalized to more complex situations like those in which noise or external 
forces affects the behavior of the system. 
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